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In this work we applied a quantum circuit treatment to describe the nuclear spin relax- 
ation. From the Redfield theory, we were able to describe the quadrupolar relaxation 
as a computational process in the case of spin 3/2 systems, through a model in which 
the environment is comprised by five qubits and three different quantum noise channels. 
The interaction between the environment and the spin 3/2 nuclei is then described by 
a quantum circuit fully compatible with the Redfield theory of relaxation. Theoretical 
predictions are compared to experimental data, a short review of quantum channels and 
relaxation in NMR qubits is also present. 
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1 Introduction 

The analysis of physical systems in terms of information processing [TJ [2| can be very pro- 
ductive and complementary to the conventional description of nature in terms of forces and 
energy. A good example is the famous Maxwell's demon paradox, only solved using arguments 
from information theory [3]. Recent works have applied this approach to gain novel insights 
into cosmology and fundamental physics [U [2j |4l [5l [3] . In biophysics, Engel et al. [6] found 
evidence that photosynthetic plants employ a kind of quantum search algorithm to efficiently 
capture the energy of the Sun. Furthermore, the view of physical phenomena as quantum 
computational processes can be sometimes more adequate to design algorithms to simulate 
physical systems. One open problem in quantum science is the the problem of simulating 
open systems on quantum computers, considering relaxation phenomena in the context of 
information processing could yield some insight into this problem. 

Nuclear Magnetic Resonance (NMR) is a well established technique used in Physics, Chem- 
istry, Medicine, and Biology. In the last decade, it has also been used as an experimental 
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method for many Quantum Information Processing (QIP) implementations [71 18] and in the 
study of fundamental aspects of quantum mechanics [9l[T0]. While N coupled spin 1/2 nuclei, 
either in solid or liquid state, have been extensively used to process information of N qubit 
systems, quadrupolar nuclei with spin / > 1/2, which account for about three quarters of nat- 
urally occurring magnetic nuclei, have also been used to process information of log2(2/ + 1) 
equivalent qubit systems in liquid crystals [HI [121 (El El El E] , solids [17] and more recently 
in GaAs nanodevices [HI [H [20l|2T] . 

Here we have used the quantum information view to describe the relaxation of NMR 
qubits. In this work, we describe a NMR relaxation process using the "language" of quan- 
tum circuits, following the approach of quantum information processing. For this study, the 
relaxation of a quadrupolar nucleus, originated by local electric field fluctuations, was fully 
characterized by means of basic quantum logic gates. The gates are related to some relax- 
ation parameters, such as the spectral densities and the quadrupolar coupling strength, which 
contain all the information about the relaxation process. In the present paper we demostrate 
that the relaxation of spin 3/2 nuclei can be described by a model in which the environment is 
comprised by five qubits and three different quantum noise channels. The interaction between 
the environment and the nuclei is then described by a quantum circuit fully compatible with 
the Redfield theory of relaxation. 

This paper is organized as follows: In sections [5] and [31 we present a short review of 
quantum channels and a brief description of the nuclear spin quadrupole relaxation. Sees. 2] 
and [5] contain the description of the quadrupolar relaxation as a computational process and 
a comparison between the theoretical model and experimental results. Finally, in the last 
section, some conclusion are drawn. 

2 Quantum Channels 

The dynamics of a quantum system (S) that interacts with an environment {£) can be de- 
scribed by the following Hamiltonian: 

= Hs + Ti-e + 'Hint, (1) 

where Us and He are the system and environment internal Hamiltonians, respectively, and 
Jiint is the coupling Hamiltonian between them. The set system-environment is a closed 
system and must evolve unitarily in time, according to the laws of quantum mechanics. This 
process can be represented using a quantum circuit, as may be seen in Figure [H The unitary 
evolution U = e"*^*/'' can be interpreted as a quantum computational process and therefore 
be broken into smaller quantum logical gates. Due to the action of U, after some time t, both 
S and £ become entangled and an initially pure state of S can be turned into a mixed state 

m- 

In real situations, the actual form of U , i.e. the quantum algorithm computed, depends 
on the specific processes that can take place between S and £. Although the environment 
degrees of freedom are usually very large, it is possible to model the environment with a finite 
number of qubits |22j . In fact, at most (P qubits are necessary to model the environment, 
where d is the Hilbert space dimension of S. 

Denoting p and p^nv as the initial states of S and £, respectively, one can show, by tracing 
over the degrees of freedom of the environment, that the effective evolution of p is not unitary, 
and is given by: 
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Fig. 1. Circuit model for the system-environment interaction, p and p, 
the initial state of the system and the environment. 
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represent, respectively, 
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where the set of operators Ek are the so-called Krauss operators and the condition J2k EkpE^ = 
/ must be satisfied in order to preserve the trace of the density matrix. The interpretation 
of this expression is that p is transformed into EkpEj, with probability Tr{EkpEl). The 
expression ^ is only valid if the set system-environment is initially in a non-entangled state. 
Usually, it is not the case since the constant interaction between them always produces quan- 
tum correlations. However, these correlations are destroyed upon the preparation of the initial 
state. 

Regarding the main system S as being composed by N subsystems, one can recognize two 
types of processes: The first, called Global Channel, is a process in which all the subsystems 
interact with the same environment. In this case, the operation [/ is a non-separable matrix, 
the channel can create entanglement between S and £ and, in principle, either can destroy or 
create entanglement among the subsystems. In contrast, Local Channels are those processes 
where each subsystem interacts with its own environment. In this case, entanglement among 
subsystems cannot be created since U = Ui (E) U2 <E} ■ ■ ■ <E} Un and the Eq. ([2]) has the form: 



There are many types of quantum channels, for example [22], the generalized amplitude 

damping (GAD) channel, the phase damping (PD) channel, the bit-flip, phase-fiip and depo- 
larizing channels. We will briefly review now two channels - GAD and PD - which are useful 
to describe NMR relaxation. 

2.1 Generalized Amplitude Damping 

The generalized amplitude damping channel describes dissipative interactions between the 
system and its environment at finite temperatures |22| . It can be decomposed into two 
processes. The first, denoted here as Ai^o, is a process which a qubit in the excited state |1) 
decays to its fundamental state |0) with some probability 7. In the second process (^o-s-i), 
the qubit is excited from the ground state with probability (1 — 7). The former process is 
called amplitude damping and describes a qubit in contact with a reservoir at temperature 
T = K. To introduce finite temperatures, one have to consider that the process Ao-^i occurs 
with probability V and .4i_s.o occurs with probability 1 — V, where V is the probability of 
finding the system, at thermal equilibrium, in its ground state. 
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The Krauss operators for the GAD channel are described in Eq. (l4])-([7]) [22]. The circuit 
models for the processes .Aq^i and ^i^o are shown in Figures ([2]) and ([3]). The circuit model 
for GAD can be constructed combining both processes (see Figure ([4])). 
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Fig. 2. Quantum circuit description of the process Ai^t^. The initial state of the system is 
represented by p while the environment is initialized in the pure state |0). The notation Ry{a) 
represents a a = 2 arcsin(.y/7) rotation around the axis y. 
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Fig. 3. Quantum circuit description of the process Ao^i- The initial state of the system is 
represented by p while the environment is initialized in the pure state |0}. The notation Ry{P) 
represents a /9 = 2 arcsin(.y/7) + n rotation around the axis y. 
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Fig. 4. Quantum circuit description of the generalized amplitude damping channel. The initial 
state of the system is represented by p while the environment is initialized in the pure state 
|0> (g) l*}^,, where |<I')eg = VV\0) + 0"^|1>. 
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2.2 Phase Damping 

The phase damping channel describes the loss of coherence without loss of energy [22] • In 
this channel, the relative phase between |0) and |1) remains unchanged with some probability 
A or is inverted {4> ^ (f) + ir) with probability 1 — A. The states of the computational basis 
do not change under this process. However, superpositions in the computational basis can 
get entangled with the environment. Thus, this channel does not change the probability of 
finding the qubit in the state |0) or |1), but it destroys all coherences between them. The 
quantum circuit model for this channel is shown in Figure ([5]) and their Krauss operators are 
given by [22] : 
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Fig. 5. Quantum circuit description of the phase damping channel. The initial state of the system 
is represented by p while the environment is initialized in the pure state |0). The notation Ry(9) 
represents a 6 = 2arccos(2A — 1) rotation around the axis y. 



2.3 Block Equation and Quantum Channels 

In general, two different processes occur simultaneously during the relaxation in NMR systems: 
the transverse relaxation and the longitudinal relaxation. The first leads the disappearance 
of the nuclear magnetization components {M^ and My) that are perpendicular to the strong 
static field applied along the z direction. This process causes decoherence without energy 
exchange with the environment. The longitudinal relaxation leads the ensemble of nuclear 
spins to return to its equilibrium state. This process is related to transitions between the 
nuclear Zeeman energy levels. Unlike the transverse relaxation, the longitudinal relaxation is 
a mechanism where the system exchanges energy with the environment. In order to show an 
example of how the longitudinal and transverse relaxations can be related to the quantum 
channels, let's consider the relaxation of a single nuclear spin 1/2 in the initial state 
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where Mj* is the initial magnetization along j 
channels on ((TU| leads to: 



X, y, z axis. The application of GAD and PD 
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One can show that (fTT|) - (fT3|) reproduce the solution of phenomenological Bloch equations 
[7] if 7 = 1 - e-*/^i, A = (1 + e-*/")/2 and a = 2TiT2/{2Ti - T2), where Ti and T2 
are, respectively, the longitudinal and transverse relaxation times [7]. Consequently, one can 
combine the circuits (|4]) and (O to design a quantum circuit that describes the relaxation of 
a single spin 1/2 and is completely equivalent to the Bloch equations description. The circuit 
we have derived is shown in Figure (j6|). 
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Fig. 6. A quantum circuit model equivalent to the Bloch equations description of relaxation for a 
single spin 1/2. The circuit can be constructed combining figures l(4]l and JSj. 



3 Nuclear Electric Quadrupole RelsLxation 

Nuclei with spin / > 1/2 have an asymmetric charge distribution and thus possess an electric 
quadrupole moment Q. Therefore, a quadrupolar nucleus can interact with either magnetic 
fields or electric field gradients [331 121] ■ Usually in liquid state, rapid molecular tumbling 
tends to completely average the electric field gradients and consequently, this interaction is 
not observed in the spectrum. However, even when any quadrupolar splintting is observed, 
the relaxation due to electric field gradient fiuctuations can significantly contribute to the 
longitudinal and transverse relaxation. In many cases of interest, the averaging is not complete 
and an residual axially symmetric electric field gradient can be effective. In this cases, the 
dynamics of a quadrupolar nucleus is described by the Hamiltonian [7]: 
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where and ujq are, respectively, the Larmor and quadrupole frequencies. For a spin 3/2 
system, such as the ^^Na nuclei, this Hamiltonian gives rise to four unequally spaced energy 
levels, originating an NMR spectrum containing three lines, corresponding to transitions 
between adjacent levels. The energy states |3/2), |l/2), | — 1/2), and |— 3/2), can be labeled 
as |00), |01), |10); and |11), which correspond to a two-qubit system. 

Nuclear relaxation is caused by the interaction between the nuclear spin and random 
electromagnetic fields generated by the environment. The relaxation can occur through several 
mechanisms, depending on the type of interactions existing in the spin system [53]. When the 
quadrupolar coupling strength is much higher than the dipolar fields, the loss of coherence 
and energy dissipation can then be considered to be exclusively due to electric field gradient 
fiuctuations [231 1241 125] . Under pure quadrupolar relaxation mechanism, one can show, using 
the Redficld formalism ([23l[26]), that the time evolution of each element of a spin 3/2 density 
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matrix is given by [27] the equations (fT5|) -(|24 |) . The superscript eq denotes the thermal 
equilibrium state and the index values i, j = 0,1,2,3 corresponds to the states |3/2), |l/2), 
| — 1/2), and |— 3/2), respectively. 

Poi [t) = ^[Poi (to) + P23 (to) 

+ (POI (to) - P23 (to)) e-2C7.(t-to)]g-C(Jo + J.)(*-*c)^ (I5) 

P23 (t) = ^[Poi (to) + P23 (to) 

- (POI (to) - P23 (to)) e^2C:j.(t-to)],-C(,/o + .A)(*-*o)^ (16) 

P02 (t) = -[P02 (to) + P13 (to) 

+ (P02 (to) - P13 (to)) e-2C^i(*-*°)]e-^(-^«+-^^)(*-*''', (17) 

Pl3 (t) 2 [^02 (to) + Pl3 (to) 

- (P02 (to) - P13 (to)) e-2C:j.(t-to)],-C(,/o + ./.)(*-*o)^ (18) 

Pi2(t)-Pi2(to)e-'='('^^+'^^)(*-*''\ (19) 
Po3(t)=Po3(to)e-^(-'^+-^^)(*-*''', (20) 

Poo yt) ~ Poo - 4^^!"^ 



where 



-ROe-^CMt^to) _ ^0g-2CJi(t-to)]^ (21) 
Pii(t)=Pil + ^[^?e-2C(./i+./.)(t-to) 
+i?°e-2CJ2(t-to) _ ^0g-2CJi(t-to)]^ (22) 

P22 (t)=pl + i[i??e-2^('^^+'^^)(*-*°) 

_^0g-2CJ.(t-to) + ^0g-2CJi(t-to)]^ (23) 

P33(t)=P3l-i[i??e-2^('^^+'^=)(*-*°' 
+i?°e-2^-^^(*-*°) + i?^e-2^^i(*-*°)]. (24) 
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(25) 
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The parameters Jq, Ji and J2 are the spectral densities [2 8) and C is a proportionality 
coefficient that can be determined using the quadrupolar coupling strength [29j . 
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4 Quantum Computing Model for the Quadrupolar Relaxation 



In the work of Auccaise et al. [27], it was shown that the longitudinal relaxation of each 
qubit in a spin 3/2 nucleus under pure quadrupolar relaxation are related to the spectral 
densities J2 and Ji, respectively, and occurs independently. These features suggest that the 
longitudinal relaxation in this case can be described by two dissipative GAD (see section [O]) 
channels acting on each qubit separately. 

In contrast, the loss of coherence depends on the three spectral densities. By inspection 
of Equations ()15p - (j20p . one can sec that only the coherence elements poi, po2 , P13 and P23 
depend on Jq, suggesting that the phase damping channel in this case does not destroy the 
coherence between the states |00) and |11) and between |01) and |10), in other words we 
can say that the Bell's states form a free decoherence subspace for such channel. A possible 
quantum channel to reproduce this feature is a channel where the relative phase between 
|0) and |1) in both qubits remains unchanged with some probability A or are simultaneously 
inverted with probability 1 — A. We will call this channel as Global Phasing Damping (GPD). 
The quantum circuit model for the GPD channel is presented in Figure ([7]) and its Krauss 
operators arc given by: 
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Therefore, the quadrupolar relaxation process can be described by a model with three 
different channels: two local channels related to dissipation of energy and one global channel 
related to loss of coherence without energy exchange. From the above statement, we can 
conclude that the quadrupolar relaxation can lead to sudden death of entanglement. Since 
the Bell states are not affected by the phase damping channel, the loss of coherence in such 
states occurs due to two independent GAD channels. As demonstrated in [50], under this 
situation all states undergo sudden death of entanglement. One can show, using Eq. ^ and 
Q, that under the action of these channels, the evolution of each density matrix clement is 
given by: 
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poi(i) - [{l-lA{l-VA))poi{to) 



+ hAVA)p23ito)W'^ - lBi2\ - 1), (30) 

P23{t) = [(7A(l-7'A))poi(to) 

+ (1 - lArA)p23ito)Wl - 1b{2X - 1), (31) 

P02(i) = [(l-7s(l-7'B))po2(io) 



-(7B^B)pi3(to)]v/T^(2A - 1), (32) 



Pl3{t) = [7b(1 - '^B)/302(io) 

+ (1 - lBVB)pi3{to)W'^ - 7a(2A - 1), 
Pl2{t) = Pi2(io) Vl - 7A - IB, (33) 



/003(i) = Po3(io) ^ 7A ^ 75, (34) 

Poo(i) = (1-7A(1-7'A))(1-7B(1-PB))/Ooo(io) + 
(1 - 7a(1 - 7'A))(7B7'B)/Oii(to) + 
(7a7'a)(1 - 7s(l - ^B))/522(to) + 

ilAVA)hBVB)p33{to), (35) 
/9ll(i) = {l-jA{l-VA)){jB{l-VB))poo{to) + 

(1 - 7a(1 - Pa))(1 - 7B7'i3)pii(to) + 

(7A^A)(7B(l-7'B))p22(to) + 

ilArA)il-lBVB)p33ito), (36) 

P22(i) = (7A(l-7'A))(l-7B(l-7'B))poo(to) + 
{-fAil ~ VA))hBVB)pil{to) + 

(1 - 7a7'a))(1 - 7s(l - rB))p22{to) + 

{l-lAVA){lBVB)p33{to), (37) 
P33(i) = {lA{l-VA)){lB{l~VB))poo{to) + 

(7a(1 - Va)){1 - lBVB)pii{ta) + 

{l~lAVA){lB{l-VB))p22{h)+ (38) 
(1 - 7a7^a)(1 - lB'PB)P33{to)- 



Comparing ((TS|) - (P^ and (pn|) - (|38p . one can find the conditions for which the combined 
action of the three quantum channels is equivalent to the Redfield theory for the quadrupolar 
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relaxation: 7^ = l-e-^^-'''^ 75 = l-e-^'^-'^\ A = ^ {l + e-^^°*) and T'a = "Pb = 1/2. Note 
that the last condition implies that the equilibrium state is a maximally mixed state, which 
corresponds to the equilibrium state at infinite temperature. Nuclear spins under Zeeman 
interaction at room temperature correspond to a low polarized system, ss 1/2 for all 
spins). Thus this requirement is quite reasonable, as we will show comparing the theoretical 
prediction to the experimental data in the next section. Furthermore, it is important to 
emphasize that a high temperature approximation is implied in the solutions of Rcdfield 
equations (|T5|) -(|24 |) . 

Therefore, a quantum circuit model for the quadrupolar relaxation can be constructed 
combining the quantum circuits ([4]) and ([7]). The circuit equivalent to the relaxation is shown 
in Figure ([9]) and involves 7 qubits. The first two qubits, from top to bottom, correspond 
to the qubits of the quadrupolar spin. The remaining five qubits correspond to the environ- 
ment, which is initialized in the state |000) ® 1*1') eg, where \^)eq 



(v^b|0) + \Jl — Vb\^))- The relations between the angles 
densities are summarized in the table ([T|). 



{^/V2\()) + VT^V^\l))® 
ctA.B, I3a,b and the spectral 



Table 1. Relation between 8, ctA,B^ Pa,b and the spectral densities. 



a A = 2 arcsin ( \/l — e^2CJ2t 
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/3a = 2 arccos ( Vl - e^^cj; 



13b = 2 arccos (^VJ^e^^cJlt 
e = 2 arccos f2e-<^-^o*) 
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Fig. 7. Quantum circuit description of the global phase damping channel. The first two qubits 
(on the top of figure) correspond to the qubits of the quadrupolar spin. The environment is 
comprised by only one qubit initialized in the pure state |0). The notation Ry(8) represents a 
8 = 2 arccos(2A — 1) rotation around the axis y. 
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5 Experiment 

The experiments were carried out using a 9.4 T - VARIAN INOVA spectrometer on ^^Na 
nuclei in a sample of a lyotropic liquid crystal system (Sodium Dodecyl Sulfate = SDS) 
|31j . Several pseudo-pure states were prepared using numerically optimized pulses known 
as Strongly Modulated Pulses (SMP)[32l [33] and the evolution of each density matrix ele- 
ment was monitored by quantum state tomography [T2j El [34]. More specifically, it was 
experimentally followed the time evolution of the pseudo-pure states [7] : 

Ppp. = ^^5 + f|V')(^l, (39) 

where e ^ 10^^ is the nuclear polarization at room temperature and Ii/j) corresponds to |00), 
|01), |10), 111) i(|00) + |01) + |10) + 111)), and ^(|00) + |11)). 

The basic experimental scheme presented in Figure ([8]) consists of: a state preparation 
period performed using SMP technique |321 133j ; a variable evolution period where relaxation 
processes take place; a hard RF pulse with the correct phase cycling and duration to execute 
quantum state tomography via coherence selection |34| . For off-diagonal elements, a tt pulse 
was added in the middle of the evolution period to refocus the Bq field inhomogeneities. 
Since the quadrupolar evolution is not refocused by the tt pulse, evolution periods multiple 
of 2-K /loq were used. The spectral densities were experimentally determined in reference |27| 
as being Jq = (14 ± 1) x IQ-^s, Ji = (3.4 ± 0.4) x lO'^s, and J2 = (3.7 ± 0.3) x lO'^s and 
C = (1.2 ±0.1) X lOiOs-2. 

Using the experimental spectral densities and the initial density matrix as inputs, we have 
simulated the evolution of density matrix (|39l) using the 7-qubit circuit model shown in ([9]) 
and then compared our simulation with experimental data. The comparisons for populations 
of deviation density matrix (Ap = ppps — //4 [7|) corresponding to the states |01), |10), |11), 
^(|00) + |11)) and i(|00) -I- |01) + |10) + |11)) arc shown in Figure [TO] The results for the 
non-diagonal elements corresponding to an uniform superposition are presented in Figure 
pip . In all cases, the experimental behavior could be well reproduced using the circuit model 
([9]) . The good agreement between the theoretical and the experimental results shows that the 
electric quadrupolar nuclear relaxation can be viewed as three different channels and that the 
circuit derived in this work provides a good quantum computation model for the relaxation. 
Small discrepancies between the experimental and theoretical results arc due to experimental 
errors in the determination of the spectral densities and initial density matrices. 

6 Conclusions 

In this work, we have studied the relaxation of a nuclear quadrupolar system, originated by 
local electric field gradient fiuctuations, using the language of quantum circuits and following 
the approach of quantum information processing. Wc were able to describe the quadrupolar 
relaxation by a model in which the environment is comprised by five qubits and three different 
quantum noise channels acting on the quadrupolar system. The interaction between the 
environment and a spin 3/2 was described by a quantum circuit fully compatible with the 
Redfield theory of relaxation. The theoretical predictions were compared to experimental 
results. The good agreement between the theoretical and the experimental results shows 
that the model derived in this works provides a good quantum computation model for the 
relaxation. 
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Fig. 8. Scheme of the pulse sequence used for probing the relaxation of the individual density 
matrix elements. The initial states are prepared with the SMP technique |32) . The state relaxation 
takes place during a variable evolution period and, finally, a hard RF pulse with the correct phase 
cycling and duration is applied to execute quantum state tomography via coherence selection |34) . 
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Fig. 9. Quantum circuit description of the nuclear spin 3/2 electric quadrupole relaxation. The 
first two qubits (on the top of figure) correspond to the qubits of the quadrupolar spin. The 
remaining five qubits corresponds to the environment which is initialized in the state |000) C>5 |"I')eg 
where |$>e, = {VPI\0) + V^-T^aW) ® -(v^lO) + V^-Vb\1))- 
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Although many studies concerning open systems have been carried out for decades, just 
few studies consider relaxation phenomena in the context of information processing. We 
believe that further efforts on this direction could bring a better understanding about open 
systems as well to yield some insight into the problem of simulating open systems on quantum 
computers. 
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Fig. 10. Combined experimental data ( -poOi ▼-Piii B-P22i ^ -P33 ) and the corresponding theoret- 
ical prediction (solid hnes) obtained from the circuit ^ for various input states. The experimental 
data were obtained from previous experiments |27| . 
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Fig. 11. Combined experimental data (symbols) and the corresponding theoretical prediction 
(solid lines) obtained from the circuit ||9ll and starting from the state ^(|00) + |01) + |10) + |11)). 
The experimental data were obtained from previous experiments |27| . 
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